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We present a theoretical analysis of the phase diagram of two-component bosons on an optical 
lattice. A new formalism is developed which treats the effective spin interactions in the Mott and 
superfluid phases on the same footing. Using the new approach we chart the phase boundaries 
of the broken spin symmetry states up to the Mott to superfluid transition and beyond. Near the 
transition point, the magnitude of spin exchange can be very large, which facilitates the experimental 
realization of spin-ordered states. We find that spin and quantum fluctuations have a dramatic effect 
on the transition making it first order in extended regions of the phase diagram. For Mott states 
with even occupation we find that the competition between effective Heisenberg exchange and spin- 
dependent on-site interaction leads to an additional phase transition from a Mott insulator with no 
broken symmetries into a spin-ordered insulator. 

PACS numbers: PACS 



I. INTRODUCTION 

Recent observations of the superfluid to Mott insulator transition in a system of ultracold atoms in an 
optical lattice open fascinating prospects for studying many-body phenomena associated with strongly 
correlated systems in a highly controllable environment PjMH- ^'^'^ instance, theoretical studies have 
shown that, with spinor bosonic or fermionic atoms in optical lattices, it may be possible to observe 
complex quantum phase transitions to realize novel superfluidity mechanisms 0, and to probe one- 
dimensional systems exhibiting spin charge separation Q. 

Recently, Duan et al. proposed a technique to implement interacting spin-i Hamiltonians using 
ultra-cold atoms, opening the door to controlled studies of quantum magnetism. In this approach the 
two-state bosonic or fermionic atoms are confined in an optical lattice where spin-dependent interactions 
and hopping are controlled by adjusting the intensity, frequency, and polarization of the trapping light. 
Deep in the Mott phase the motional degrees of freedom are frozen out and the remaining spin degrees 
of freedom are coupled by an effective Heisenberg exchange. In refs. an effective spin hamiltonian 

was derived by perturbation theory for the case of a single atom per site and the limit of small tunneling. 
However, in practice Mott states with more than one atom per site are also of considerable interest and 
may exhibit richer phase diagrams. Furthermore, spin effects are expected to be important, and even 
stronger, at larger values of the tunnelling, where perturbation theory fails. For example, an important 
question that cannot be addressed by the perturbative treatments is how spin affects the transition into 
a superfluid phase and the properties of the superfluid phase itself. 

In this paper we first extend the earlier approaches to the case of Mott states with general integer 
occupation. We find that at even fillings the competition between on-site interactions and nearest neighbor 
spin exchange leads to a transition from a spin ordered Mott state to one with no broken symmetries. 
Then we present a theoretical framework, which is non-perturbative in the tunneling and allows to 
describe both the superfluid and insulating phases in two-component systems. Using this approach we 
determine the phase diagram for a density of one atom per site. We find that the spin-ordered states 
persist up to the superfluid transition. In this region the critical temperature for spin ordering can be 
large, facilitating experimental realization of these phases. The z-antiferromagnetic state in particular, 
enjoys a negative zero-point energy which extends its domain far beyond the mean field prediction for 
the Mott phase. The transition between this state and the superfiuid is found to be first order in contrast 
with the standard superfluid-insulator transition. 

Before proceeding, we note that spin Hamiltonians can also be simulated by controlled collisions via fre- 
quent time-dependent shifts of the lattice potentials @. Compared with that method, the spin-dependent 
tunneling may have certain experimental advantages since it implements the desired Hamiltonian directly, 
thus circumventing imperfections and errors associated with rapid perturbations due to the lattice shifts. 
We also note the recent studies on quantum magnetism induced via magnetic dipole interactions of the 
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condensed atoms [lOj. The present approach results in much larger interaction strength per atom, and 
also allows for more flexible control over interaction properties. 

The paper is organized as follows. In section |nl we describe the Hubbard model for two bosonic species 
on an optical lattice, which serves as our starting point. In section IHII the perturbative approach of refs. 
0, 111 is extended to arbitrary integer filling and the insulating phase diagram is investigated. In section llVl 
we present the mean-field description of the SF-MI transition in two component systems. The analytical 
predictions of a variational approach are compared to the results of a numerical mean- field analysis. In 
section a theoretical framework is developed that incorporates the effect of quantum fluctuations and 
treats the magnetic interactions in the Mott and superfluid phases on an equal footing. In Section IVII 
this framework is used to analyze the full phase diagram for one atom per lattice site. The relevance of 
the present results in the light of realistic experiments is discussed in the concluding paragraphs of the 
paper. 



II. THE MODEL 



We consider a system with two species of atoms or equivalently, atoms with two relevant internal states. 
The two species shall be denoted by the second quantized bosonic operators a and b. We assume that the 
two species are trapped by independent standing wave laser beams through polarization (or frequency) 
selection. Each laser beam creates a periodic potential in a certain direction Vaasiv^ikar), where is 
the wavevector of the light and cr = a, 6 is the species index. Throughout this work, we assume that the 
laser beams are orthogonal, creating either a square lattice in two dimensions or a cubic lattice in three 
dimensions. For sufhciently strong periodic potential and low temperatures the atoms will be confined to 
the lowest Bloch band. The low-energy Hamiltonian is then given by the Bose-Hubbard model for two 
boson species: 

H ^ - Y,ta {a\a, + H.c) - {blbj + H.c) 

(ij) iv) 
i 

+ ^ ^ VaUatiUai - I) -'^Hallai (1) 
ia—a.b ia 

Here {i, j) denotes the near neighbor sites, a^, bi are bosonic annihilation operators respectively for bosonic 
atoms of different spin states localized on site i, Uia = a} ai^, Uib = foj For the cubic lattice, 

using a harmonic approximation around the minima of the potential Q, the spin-dependent tunneling 
energies and the on-site interaction energies are given by ta(b) ~ (i'^/4) Va(b) exp[— (7r^/4) (wa(&)/£'fl)^^^], 
U ~ (8/7r)^/^(fcaa^fc)(i?flvj^j)^/''. Here Va,b is the depth of the optical potential for species a and 6, 
Vab — ^VaVb/{vl^^ + "^fc^^)^ is the spin average potential in each direction, Er = h?'k'^ /2m is the atomic 
recoil energy, and aab is the scattering length between the atoms of different spins. The intra-species 

interaction is given by Va(b) ~ (S/tt)^/^ (fcaa(;,)) f^'flw^^f,)] {a,a{b) are the corresponding scattering 

lengths). Furthermore, the magnitude of the interspecies interaction U can be additionally controlled by 
shifting the two lattices away from each other, which opens a wide range of U /Va to exploration. Note 
that spin-dependent tunnelling can be easily introduced by varying the potential depth Va and Vb with 
control of the intensity of the trapping laser. We should also point out that the two atomic states generally 
have different energies (/Iq ^ /ib in0. In the spin language, this translates to a magnetic field in the z 
direction. However since there is essentially no transfer between the two populations, the experiment is 
performed with fixed magnetization and the chemical potentials can be set to fix this magnetization. 

In this paper, we address primarily the case in which the total filling is commensurate with the lattice 
and the two species have equal density. A transition from a superfluid to a Mott insulator is expected, 
as in the usual case of a single species. However, in this system, magnetic order, associated with the 
pseudospin degrees of freedom (boson components), may occur as well. 
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III. DEEP MOTT PHASE: EFFECTIVE SPIN HAMILTONIAN 



To illustrate the magnetic orders that can arise it is instructive to begin deep in the Mott insulator 
in the limit ta^b << U,Va^b: where the hamiltonian ^ can be simplified considerably. The low energy 
Hilbert space in this case contains states with a particular integer occupation on every site. However 
there is a remaining degeneracy associated with the spin (boson component) degrees of freedom. The 
degeneracy can be removed by an effective hamiltonian acting within the low energy subspace. This was 
done previously in Refs 7 , S] for the case of a single atom per site, using second order perturbation theory 
in the hopping parameters. The result is 

iJe// ^J.J2 S^^l - H^StS^ + S^S^) - hY^Sl (2) 

{ii) (ij) i 

Here 1 1 ) and | J, ) represent sites occupied by the a and b atoms respectively and the couplings are given 
by: 

J ^ /b+tl ^tl Atl 



J± 



U Va Vb 

4:tatb 
U 

2tl 2tl 



We assume that the induced ordering field h can be cancelled by an externally applied field hext- In this 
case the model obviously exhibits a transition between a. x ~ y ferromagnet for J± > Jz > to an Ising 
antiferromagnet with z-Neel order (Fig. 




FIG. 1: Schematic phase diagram of the effective spin hamiltonian ((2J, valid deep in the Mott phase. 



We now extend the discussion to the case of any integer filling of N atoms per site. To see how things 
may become qualitatively different from the singly occupied case, consider first a Mott state with two 
atoms per site. The low-energy Hilbert space of a lattice site consists of the three states 

-i=(at)2|0), at6t|0), ^(&^)'|0). (4) 

If Va.b >> U, the state a^fet |0) has much lower energy than the other two. This implies a simple 
Mott state J^^ a|6j | 0) which, unlike the Mott states in Fig. does not break any symmetries. On the 
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other hand, when Vafi is of the same order as U, all three states should be taken into account and more 
interesting phases may be possible. Therefore in the general case of N atoms per site we consider the 
regime ta,b, \Va,b - U\ « U, K.fa and K = Vf,. 

The low-energy Hilbert space of a lattice site with N atoms per site can be constructed in a similar 
way. It contains the + 1 states: 

\S,m) ^^^J== ^ |0) (5) 

where S = N/2 and m — —S, . . . , S. Obviously and act as Schwinger bosons creating a multiplet 
of pseudospin S. The spin magnitude depends on the site occupancy. It is integer for even N and half 
integer for odd N. 

Now the effective hamiltonian within the spin S subspace can be derived by second order perturbation 
theory in a straightforward generalization of Ref . 7] . The result is 

Heff = - E ["^j- (^^^j + ^'^J) + J^^^^^] +u^{s!r-hY^ s! (6) 

{ij) i i 



where the interactions are given by: 



J 



u - Va-U = Vb-U 



J. = 2^ 



U 

tl + tl 



u 



h = z{2S+iy-^j^ + K,t (7) 

Note that if we take the parameters are identical to those of the effective spin-i hamiltonian jSJl in 

the case U ^ Va — Vf,. Note also that additional terms such as (SfSj)'^ do not arise. The (5*^)^ term is of 
course just a constant for S = and therefore has no effect in this case. However, it plays an important 
role at larger value of the spin, namely for occupations > 1. 

When u can be neglected relative to or J± , the remaining terms in ^ form an anisotropic Heisenberg 
model. A standard, coherent state mean field theory is then possible, which yields x — y ferromagnetic 
order. A large positive u, in the absence of an ordering field {h = 0), acts to reduce the Sz component 
of the spins. At some point, the classical coherent states that represent fully polarized spins, become 
unsuitable descriptions of the system. In particular, at large enough u all spins will be essentially confined 
to their lowest possible states. 

When the spin is half integer (odd filling), there are two active states, at large u, corresponding 
to = ±i. The hamiltonian © then reduces to a spin-i model, but the spin interactions remain 
practically identical. Thus, the essential physics is unchanged. We expect ferromagnetic spin order, as 
for small w, only with a reduced effective moment. 

For integer spin (even filling) we expect qualitatively different behavior. At large enough u only 
the ~ state will be important. We then expect that the system is well described by | ^ ) = 
Hi \ S,m = 0) a Mott state with no broken symmetries. The transition to this state at large u from 
the X — y ferromagnet at small u is formally identical to the transition from a superfluid to a Mott phase 
in the Hubbard model of single component bosons. A direct correspondence exists between the boson 
number in the Hubbard model and in The Mott state of bosons is characterized by vanishing 
particle number fluctuations on a site. Similarly the transition here is into a state with well defined 
on each site. 

To describe the transition we note that only the three states with lowest play an important role in 
its vicinity. We therefore write a homogenous mean field ansatz: 



n [ cos(0/2) I 5, ) + e"' sin(0/2) (e'^ cos(x/2) | 5, 1 ) + e''^ sin(x/2) | 5, -1 ) ) ] . (8) 
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The variational energy in this state is given by: 

J z J z 6 

E = —S{S + 1) sin^ 9{1 + sinxcos2ry) ^■^^'^^o ^'-'^^ ^ '^sin^ - — /isin^ - cosx (9) 

8 2 2 2 2 

and we see that the minimum occurs for rj = 0,7r. The x — y order parameter is (S*^) oc sinfl = tp. 
Therefore, to find the transition to a Mott insulator we expand the energy up to quadratic order in ip 
and minimize it with respect to x- Note that the quartic term is always positive since Jz < J±. We then 
obtain the critical value of Jj_ as a function of h: 

u zS(S + 1) 

For magnetic fields h > u, the description in terms of the states {| — 1), |0), |1)} breaks down. 
Instead, a similar scheme can be carried out, using the states {|0), |2)}. This yields another lobe 
corresponding to a phase with well defined = 1. A schematic phase diagram is plotted in Fig[21 As the 
ordering field is increased, we obtain lobes corresponding to larger values of up to = 5*, where the 
spin is fully polarized. In practice the number of particles in each spin state is conserved independently. 
In other words the experiment is done with fixed z magnetization and h is used as a theoretical tool to 
set this magnetization in our model. Here we fix zero magnetiztion by setting h = 0. 




FIG. 2: Phase diagram of tlie Mott state with even number of atoms per site. The lobes mark Mott states with 
fixed = {ua — nb)/2. Outside of these lobes the system is also in a Mott state but with x — y ferromagnetic 
spin order. 



In summary, we found that the Mott phases of two-component bosons with even site filling are markedly 
different from those with odd filling. At odd filings, the Mott regions of the phase diagram are essentially 
the same as those found for single occupation (see Fig. These are all broken symmetry phases, either 
a, X — y ferromagnet or a z-Neel state. The Mott phases at even filling are sketched in Fig. |2] Most 
notably, another Mott-type transition occurs within the Mott phase, between a x — y ferromagnet and 
a non symmetry breaking Mott state. Related spin ordering transitions for spin 1 bosons have been 
discussed recently by Imembakov et al. 0| and by Snoek and Zhou [T3j| . 

The perturbative expansion leading to ((SJ breaks down as the transition to a superfiuid is approached 
and ta^b become comparable to U. The question arises, whether the phases predicted by the effective spin 
Hamiltonian still hold in this regime. More importantly, how do the effective spin interactions affect the 
nature of the transition to a superfiuid and the superfiuid phase itself? 
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To answer these questions we shall develop in the next two sections a theory which captures the effective 
spin interactions while also able to describe the transition to a superfLuid. 



IV. MEAN FIELD THEORY OF THE SUPERFLUID-MOTT TRANSITION 

The usual, single component, Mott transition of bosons is well described by mean field theory 
It is thus natural to start our treatment of the two component case with a mean field approach. In order 
to capture the superfluid phase we extend the regime considered in the previous section to allow for 
arbitrary ratios of ta.b/U. However we shall confine ourselves to the case of a single atom per site and to 
the limit U, ta.b << Va, Vt. Later we shall consider corrections due to finite intra-species interactions. 

In this lirnit, it is particularly advantageous to use a variational approach, which is equivalent to mean- 
field theory The idea is to assume a site factorizable wave function associated with hard core bosons, 
which in our case takes the form 



1$; 



(si 



2 " 



2 



sm —a'^+ cos ^b. 

cos — sm h cos — a' b] 

2 V 2 2 * ' 



(11) 



The enormous reduction in Hilbert space, made possible by neglecting double occupation, is what makes 
these states convenient to work with. Specifically, it is easy to calculate expectation values. In addition, we 
shall see that they facilitate a fluctuation expansion about the mean field theory. Generalization to include 
higher occupations is possible but would make the subsequent calculations much more complicated. Note 
that a more general mean-field ansatz would include complex weights, however it is easily verified that 
this would not improve the variational energy. 

In the Mott state, where each site is occupied by exactly one atom, the variational state simplifies even 
more 



MI , 



n 



iia/2 ■ Xi t 

e ^' sm — a' 



10) 



(12) 



It can be viewed as a pseudospin-i state with | ) = 1 1 ) and | ) = 1 1 ) . 

The onset of superfluidity is characterized by the development of an order parameter sin 6 ^ 0. More 
precisely, the superfluid order parameters of the two species in the state | $ ) are given by: 



1 • /3 f X-V 

- sm cos — - — 

2 V 2 

1 . n ■ f X + V 

— sm tl sm 

2 \ 2 



(13) 



Now a classical energy functional can be written, which is defined by the expectation value of in | $ ) 
Allowing for two sub-lattice order the energy function is: 



E 



zta . a - a (xa~Va\ ( Xb-Vb 

sm &A sm B cos cos 

4 \ 2 ) 

. a . r, ■ fXA + VA\ . f XB +VB 

smffA smt'R sm sm 

4 \ 2 J \ 2 

-|-^(cos 9a + cosOb), 

o 



(14) 



where z is the lattice coordination number. In the superfluid phase this function is minimized when both 
cos((xi — rii)/2) — 1 and sin((xi — Vi)/'^) — 1; which implies Xi = Vi = Tr/2. The remaining degree of 
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freedom 9 is uniform on the lattice and found by minimizing 



E{9) = +tb)sin2(6') + ^cos6' (15) 



The result is 



aC0s(-tc/(ia +4)) ta+tb>tc ^^^"^ 

where tc = U/2z. We thus find a transition to a Mott insulating state for ta + h < tc as illustrated by 
the circles in Fig. 13 This constitutes a straightforward generalization of the standard transition for a 
single species. 

By assuming the variational state pi|l , we neglected contributions from states with multiply occupied 
a or & bosons. To determine effects arising from the finite magnitude of the intra-species interaction we 
use a numerical self-consistent mean field field theory of iQ). As first proposed in (IJj, the kinetic energy 
terms in the Hamiltonian are decoupled: 



MF 



i:a—a.b 



- (aliaj) + H.c^ -hY^ {bl{bj) + H.c^ + const. (17) 

(ij) (ij) 

In the homogeneous phase this leads to a sum of identical single-site Hamiltonians 

Hmf = U iua - ^ j iub - ^ ) + ^ ^ VaUa (Ua - 1) 



2/ V 2/ 2 

OL—a,b 

-^a{a^ +H.c.)--9b{b^ + H.C.) (18) 
where the decoupling fields have to be determined self-consistcntly according to 

•^a,b^ Zta,b{a{b)). (19) 

We have solved the combined set of Eqs. H18|l and (|19|l numerically by diagonalizing Hmf within a 
finite-size Hilbert space where we allow for up to M=9 bosons per species. 

We show results in Fig. ((SJ, where it can be seen that for a small ratio U/Va,b the phase diagram is 
identical to that determined variationally. As 14,6 decrease and approach U the Mott domain shrinks. 
For Va^b < U there is an instability toward a z-ferromagnetic superfluid. Since the experiment is done at 
fixed magnetization this would lead to phase separation into domains occupied only by a or by 6 atoms. 

Note that in the Mott state where the order parameters (a) and (6) vanish, the ground state of Hmf 
has precisely one atom per site but is completely independent of the relative weights of a and b atoms. 
Similarly, the variational energy in the Mott state H12|l is a constant (— 1//4), independent of the individual 
spin orientation. Thus the simple mean field approaches are unable to resolve spin order in the Mott state. 
To obtain spin order we shall in the next section consider quantum fluctuations around the variational 
mean field solutions. 



V. EFFECT OF FLUCTUATIONS: "MAGNETIC" STATES 



The situation we encountered when attempting to treat the Mott phase with the variational states is 
similar to the basic problem of frustrated quantum magnets. The classical energy of such systems, i.e 
the expectation value of the hamiltonian in a basis of coherent spin states, often contains a macroscopic 
degeneracy (see for example the review Ref. ^3). A general mechanism that can lift the degeneracy is 
"quantum order by disorder", whereby broken symmetry configurations are selected by the zero-point 
energy due to spinwaves [l7| . 



8 




2zt/U 

a 

FIG. 3; Phase diagram obtained by the decoupling mean-field theory (I17II for U — 20 and different values of 
Va.b- Note that for finite Va,b and strong asymmetry of the hopping ta,b one of the species can be completely 
depopulated in the superfluid phase (see lower left inset). Except for this case, the Mott transition always happens 
at the same parameters for both species. 



A spinwave expansion in magnetism includes the quadratic fluctuations around coherent-state mean 
field configurations. We formulate a similar expansion in fluctuations about the mean-field states Hll|) . 
As a first step we define second quantized bosonic operators, that create the appropriate Hilbert space: 



{a| \0),bl \0),albl I 0) , I 0) } ^ {a^ \n) ,al\n) ,pl\n) ,hl\n)}, 



(20) 



where | ) is the vacuum of the new bosons and | ) is an empty site. The new operators are analogues 
of Schwinger bosons in spin systems. Like the Schwinger bosons, they obey a holonomic constraint, 
namely that their total filling on a site is one. Now, we apply an orthogonal change of basis: 





( 










\^3^ J 


V 



Xi. 



sm cos 



cos Y sm Y cos ^ cos ^ 



cos ^^ 



-sinf 



■ sin sin -y 


cos ^ 



sin Y cos Y" 
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/ Ctii \ 








Pi 









In the new basis the variational state is simply a singly occupied Fock state of the -00 boson 



(21) 



(22) 



The three remaining bosons, V'l 2 3j create orthogonal fluctuations about the variational state. 

In the constrained Hilbert space of no double occupancy by the same species, the hamiltonian may 
be written in terms of the ^ bosons. Furthermore, V'oi can be eliminated using the hard core constraint 



\ 



(23) 



so that the hamiltonian is a function of only the three fluctuation operators V'i,2.3- Assuming the fluc- 
tuations are small, we expand it to quadratic order in these operators. The exact form of the quadratic 
hamiltonian depends on the variational starting point which fixes the rotation matrix H21f) . For a two 
sub-lattice variational state, the fiuctuation hamiltonian has the general form: 



H 



flue 



Gk 

T* 



*fc - tr Tk 



(24) 
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where 

*lc = i< i^U. ^-k v-ik+J, = (v-L V'L V'L) (25) 

while J^k and are 6x6 matrices which depend on the variational parameters. Finally Hfi^c is 
diagonalized by a Bogoliubov transformation to obtain the excitation frequencies ujaw and the correction 
to the ground state energy: 



AE 



J-trJ^k + ^c^.kL (26) 



With the Bogoliubov transformation at hand it should be straightforward to calculate the average occu- 
pation of the fluctuations. For consistency of our approach we require: 

3 

Y,{i^L^o.^)«l (27) 
a=l 

Let us now focus on the Mott phase. Recall that the variational energy is independent of the individual 
spin orientations, i.e. the parameters in the state H12|l . The fluctuation hamiltonian on the other hand 
will depend on the spin configuration. Before we compare the zero point energies corresponding to 
possible spin orders let us note a few general properties of the fluctuations in this case. Since the 
bosons pI and h] which create an extra particle or hole, are unoccupied in H12|l . they constitute a 

pair of orthogonal fiuctuations. The third orthogonal fluctuation is cj)^ — cos(xi/2)a|j — sin(xi/2)a2j 
which creates a pseudospin of opposite orientation. Since the classical energy is independent of the spin 
configuration we expect that 0| will not appear in the quadratic fluctuation hamiltonian. This reflects 
the fact that a local spin flip does not cost energy. 

For a uniform state with Xi = x the hamiltonian assumes a simple form 

Hfiuc = ^Mmihk + fp{k)plpu - ^-^iplhU+Pkh-k)^ (28) 

where the couplings depend on x 

fh{k) = - {ta COS^ f + 



fp{k) = -^^ (taSin2|+if,cos2|)z7k 
g{k) = z'-fi,{ta + tb)smx (29) 



The hamiltonian is diagonalized by a standard Bogoliubov transformation: 



Pk 



coshOkCk + sinh^fcd^j, 
hk = cosh6'fcdfc + sinh^fccLj, (30) 

which yields the excitation modes 



t^l.2(k) = -\/U^ ~ 2U{ta + tb)z-/]^ + (ta + tb)"^ {zjk COS X)"^ ± {ta ~ tb)z"/]^COSX 



(31) 



In addition, there is a zero mode W3(k) = corresponding to local spin flips, which reflects the macroscopic 
degeneracy at the classical level. Higher order terms in the fluctuations take into account the corrected 
potential landscape and generate a dispersion of the spin flip mode aj3(k). Since here we are interested 
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in the zero point energy, we need not go beyond quadratic fluctuations. The quantum correction to the 
ground state energy is calculated from the prescription H26() 

AEix) - ^ E [^U' - + + + *^)' (^^k cos xY -U+{ta+ tb)z7k] - I (f + f )(32) 

k 

where we have added the last term perturbatively in ta/Va- This is justified in the regime of interest 
ta, U << Va- The minimum of AE{x) occurs for x = ""/S, which corresponds to pseudospins aligned on 
the X — y plane. Note that the dispersions of the particle and hole excitations H31|l are degenerate in this 
case. Their gap vanishes when ta + th — U/2z, which marks the transition to a superfluid in agreement 
with the variational result (fTH|l . 

To check the consistency of our fluctuation expansion the local density of fluctuations in the x — y 
ferromagnet can be calculated using the Bogoliubov transformation (|30|l : 



{Ap. + t4h^ = 4 E 2 Binh20, = ly( l-^7kfa + W[/ _ \ 



Fig. ^plots the mean square fluctuation as a function of {ta + tb)/U at a constant ratio ta/h- It can be 
verified that the mean square local fluctuation is smaller than 1/4 throughout the phase diagram. This 
constitutes a posteriori justification for our expansion which relied on the smallness of the fluctuations. 
We should comment though that the occupation of the zero mode cannot be calculated at this order. If 
the X — y state is indeed stable, interactions would generate a dispersion which would lead to a finite 
local ground state occupation. 
We now consider the canted state 

I*(X)> ^\{Us'La\ + sinh'^\ n ("sin^al+cos^&l) |0), (34) 

The angle 9 parameterizes a continuous path from the z-Neel state {9 = 0) to the x — y ferromagnet {9 — 
7r/2). Since | ^{9)) is not translationally invariant, neither will be the fluctuation hamiltonian derived 
from it. An elegant way to overcome this difficulty is to apply a unitary particle hole transformation on 
sub lattice B 

Pi ^ hi (35) 

for i €z B. In the spin language this is equivalent to a tt rotation of the spins in the B sub-lattice about 
their x axis. The rotation changes the hopping terms in the hamiltonian ^ 

alaj + h.c.^ala^j + h.c. 

blbj + h.c.-^blbl + h.c, (36) 
but it also transforms | ^'(x) ) to a translationally invariant state 

l*(x))-n(cos^al+sin^&l) |0). (37) 

Our procedure can now be carried out with the new hamiltonian and the transformed state. The fluctu- 
ation hamiltonian assumes the form 



Hfiuc = y{^iplpk + hlhi,)-^sm9{ti+t2)iplhi^ + h.c.) 



'"-^ih cos^ I + t2 sin2 '-){hihl^ + h.c.) - ^(ii sin^ cos^ l)ipW-^. + h-c) (38) 
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which can be diagonahzed by a Bogoliubov transformation. In the z-Neel state {6 
energies assume a particularly simple form 



-l,2(fc) = -Wl-(— ^ 



0) the excitation 



(39) 



Note that contrary to the x — y state, the excitations are non degenerate. The gap in wi_2(fc) vanishes on 
the lines ta^b = U/2z respectively. Thus the z-Neel state is locally stable toward formation of a superfluid 
within these boundaries. There is however a dangerous zero mode W3(k;) = which may be destabilized 
by higher order terms in the fluctuation hamiltonian. This mode corresponds to 0^ which describes spin 
fluctuations {4>l^) toward the x — y ferromagnetic state. In regions where the z-Neel state is ultimately 
stable these corrections would just generate a dispersion for W3(k). 




2z(t^+t^)/U 

FIG. 4: Mean square fluctuation in the superfluid and magnetic Mott phases. The demonstration is at a flxed 
ratio ta/ti, = 0.5. 



The quantum zero-point energy of the fluctuations in the z-Neel state is given by: 

U 



E 




(40) 



The mean local fluctuation can be calculated in the same way as before 

1 1 



-E 

2N ^ 



1 



1 



'2£t67k j2 



(41) 



It is plotted in Fig. 0] which demonstrates that fluctuations about the z-Neel state are also small. 

Before we address the full Mott domain, it is instructive to evaluate the energy corrections H32|) and 
(|40|l deep in the Mott phase, where we can compare the result with the effective spin hamiltonian 
It is also much easier to evaluate the zero point energy in this limit. For ta,tb << U we can expand the 
square roots in H32() and H40|l . then perform the momentum sums exactly, with the result 



Zjta+ti,)^ Z (tl ^ tf 
AU 

2U 



2 \ Va Vb 



(42) 
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These are identical to the mean- field energies in the effective spin hamiltonian ||2J). Thus we see that 
our fluctuation analysis about the variational states captures the essential spin interactions. In the next 
section we address the stability of the spin states over the entire parameter regime to derive a phase 
diagram. 

A fluctuation hamiltonian can be derived in a similar way for the superfluid phase where we flnd the 
three excitation modes: 



wi,2(k) = — 



ia 



k + C'(fc2) 

W3(k) - C/^t^ - 1 + (k/2)2 (43) 

with t+ = ta+ tf) and t_ = ta — tt. The zero point energy correction in the superfluid phase is evaluated 
using the prescription H26|l . A discussion of the collective modes and of the nature of the superfluid phase 
is deferred to the next section. 



VI. PHASE DIAGRAM FOR 1 ATOM PER SITE 

In this section we combine the ingredients prepared in the last sections to present a phase diagram for 
a lattice with an average occupation of one atom per site. From the variational approach we found that 
the X — y ferromagnet becomes unstable towards a superfluid state when ta + ti, > U /2z. The mean field 
phase diagram was sketched in Fig Q. However the boundaries of these phases with the z-Neel state 
remain undetermined. It is the quantum zero point energy of fiuctuations that selects ordered magnetic 
states from a degenerate variational energy. 

To analyze the stability of the spin states we need to calculate the derivatives with respect to of the 
zero point energies corresponding to these phases: 

The last term is added perturbatively in ta/Va and corrects for a large but finite intra-species interaction. 
It is easily seen that the first derivative of the modes vanishes identically at the points 9 = and 
9 = 7r/2, corresponding to the z-Neel and x — y states. Consequently these states are either minima or 
maxima of the zero point energy. The second derivative at the z-Neel state is given by: 

<P^E \ U/2 / r,(l-r,^7g) + T-a(l + r,^7k^) r,(l - r,^7^) + r,(l + r,V^) \ U 

Jg^o v/T^W yr^^ ) 2z 

where we have denoted = 2zta/U and Va = Va/U . The domain of stability of the phase is obtained 
by numerically evaluating the momentum sum in (|45|l . The resulting domain of stability is of the general 
shape illustrated in Fig. Note that the phase boundaries deep in the Mott state [ta.b «U) are linear 
and coincide with the result obtained from the effective hamiltonian (|2J). However we find in contrast 
with the effective spin hamiltonian, that even for true hard core interactions Va.fc— >oo, there is a finite 
X — y ferromagnetic domain. 

Note that the Mott z-Neel domain in Fig. Q extends beyond the mean field transition to the superfluid 
which occurs at ta + tb = U /2z. As seen in Fig. El this is due to a lower ground state energy (including 
the quantum corrrection) than the superfluid. In the remainder of this section we shall examine the 
nature of the phases and transitions in Fig. [S] 

A. Metastability and hysteresis 



It is an interesting observation that over a signiflcant parameter range, quantum fluctuations favor the 
z-Neel state even where its variational energy alone is higher than that of the superfluid. What kind 
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FIG. 5: Phase diagram including quantum fluctuations. 




0.2 0.4 0.6 0.8 1 1.2 



2z(t,+t^)/U 

FIG. 6: The energy, including quantum fluctuations, of the z-Neel state for ta/ti, = 0.1. 

of transition then, is marked by the lines tafi — U/2z, where the z-Neel state finally becomes mistable? 
It could be one of the two: (i) A first order transition into the superfluid state or (ii) a second order 
transition into a supersolid, namely a superfluid that retains Ising order. 

We shall see that the former indeed occurs, but this is not immediately obvious. Consider the excitation 
modes H39|l . Since only one of them becomes gapless on the transition lines ta,b = U /2z, one might guess 
that these lines mark the formation of a supersolid. However we now show that at the classical level the 
supersolid is unstable to formation of a uniform superfluid. 

A variational state describing a supersolid is given by 

f 9 9 \ (9 9 \ 

I ) = n «1 ( 2 ^ """^ 2 ) n ^i'^ 2 ^ 2"^ ) ^^^^ 
ieA ^ ^ ieB ^ ^ 

When 9 — TT this is just the z-Neel state. A superfluid component is added to the Neel order for 9 < n. 
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To assess stability we parameterize, a continuous deformation of this state toward the uniform superfluid: 

\m) = n 

i 

where Xi — S, tt/2 — S ior i E A, B and 6 E [0, tt/2]. The variational energy in this state, calculated using 
Eq. (Hil) is: 

E(0, J) = - ^^^^ sin 26-^ sin 2d sin S (48) 
16 16 

The second derivative {dE/d6)g^Q is negative, indicating that the supersolid is unstable. We thus es- 
tablish a first order transition between the 2-Neel state and the superfluid. We should point out that 
at higher order, quantum fluctuations can the change potential manifold. In particular, they can make 
the supersolid phase locally stable. This interesting possibility can be checked, for example by quantum 
monte carlo simulations. However, at our level of approximation there is a first order transition directly 
to the uniform superfluid. 

An important implication is the presence of a hysteresis region. Consider a change of system parameters 
from the superfluid to the z-Neel state along route A in Fig. |S1 In the region where the superfluid is 
meta-stable, it may take an excessively long time to nucleate the z-Neel state. The system is thus likely 
to remain in the superfluid state until its line of meta-stability is crossed (dashed lines in Fig[5l). Passing 
the same route in the opposite direction the z-Neel state will of course persist up to the transition line. 
Along routes B and B' matters are qualitatively different. Since the transition into the x — y Mott state 
is continuous, there is no hysteresis there. The transition from the x — y or superfluid to the z-Neel state 
is first order but without a significant hysteresis region. In fact it is well described as a transition from 
easy-axis to easy-plane anisotropy in the Heisenberg model 



sm 2 (cos y a] + sm y 6] j cos - [sm — + cos —al fe] 



(47) 



B. The superfluid and the x — y ferromagnet 

The uniform superfluid phase deserves a closer examination. In several aspects, it is different from the 
single component case. Most importantly, the two component superfluid is intimately related to the x — y 
ferromagnet. 

In terms of the bosonic operators, the x — y ferromagnetic Mott state sustains an order parameter 
(a^b) 7^ 0. In this sense it can be viewed as a counter flow superfluid It can support supercurrents 
of relative motion characterized by a gradient of the relative phase between a and b atoms. In the 
spin language this is simply a gradient of the spin orientation on the x — y plane. A magnetic field 
gradient, h(x) = h^x in Q, would twist the spin configuration at a constant rate inducing AC super- 
counterflow of frequency uj = h^L. The goldstone mode associated with this order is a spinwave, which 
describes fluctuations of the relative phase between the two components. We should note again that in the 
quadratic fluctuation hamiltonian H24|l . spinwaves are dispersionless. This is a direct consequence of the 
spin degeneracy at the classical (variational) level. However, it does not indicate truely vanishing spinwave 
stiffness. Indeed, higher order terms in the fluctuations generate a finite linear spinwave dispersion. This 
can also be understood from the effective spin hamiltonian l^l, which obviously has a finite spin stiffness 

oc Jj_. 

As system parameters are varied across the transition at ta + tf, = U/2z, one of the gapped particle- 
hole fluctuations of the x — y Mott state condenses, marking the formation of the two superfluid order 
parameters (a) and (6). As pointed out in ref. 0|, in the x — y Mott state only the relative phase is 
fixed while the average phase of the two components is disordered. In the superfluid the average phase 
orders as well. Accordingly we find two linear gapless modes in the superfluid. (|43|l . corresponding to 
an in-phase fluctuation of the two components (a;i(k)), and a relative phase, spinwave fluctuation a;2(k). 
From these considerations it is obvious that the universal aspects of the transition will be identical to the 
standard, single component Mott transition tl8l |. 
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VII. DISCUSSION AND CONCLUSIONS 



It is interesting to consider the present results in light of the current experimental possibilities. First 
of all, we have shown that magnetic phases are robust in the sense that they persist up to the transition 
into the superfluid state. Specifically, spin-ordered states appear even near the boundary of this phase 
transition. Spin-exchange interactions in this regime are most easily accessible experimentally since the 
relevant energy scales are largest and comparable to on-site interaction. In this regime magnetic phases 
are therefore relatively insensitive to perturbations due to e.g. inhomogeneous magnetic field variations. 
Secondly, we note the existence of several metastable states in this regime, indicating that the system 
is likely to display interesting dynamics as the optical potential is lowered across the transition point. 
In particular, hysteresis and abrupt changes in the state of the system can be expected. At the same 
time, our results indicate that spin-ordered states are qualitatively different for odd and even numbers of 
particles per site. Both are likely to be observable in any realistic realization, since the inhomogeneous 
trapping potential typically leads to domains with different occupation. 

Finally, it is important to note that detection of the complex states, of the type discussed in this paper, 
presents an interesting challenge in its own right. It turns out that the quantum nature of strongly 
correlated magnetic states can be revealed by spatial noise correlations in the image of the expanding 
gasT^. Specifically, atoms released from a Mott-insulating state of the optical lattice display sharp 
(Bragg) peaks in the density-density correlation function as a consequence of quantum statistics and 
such peaks can be used to probe the spin ordered Mott states proposed for two component bosons. 

In summary, we presented a theoretical analysis of the phase diagram of two component bosons on 
an optical lattice. We extended earlier treatments which were valid only deep in the Mott phase toward 
the MI-SF transition and beyond and were thus able to map a complete phase diagram. In addition we 
identified a transition into a Mott phase with no broken symmetries, which occurs only at even fillings. 
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